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I.  INTRODUCTION 


The  chemical  reaction  steps  involved  in  the  combustion  of  solid 
propellant  are  poorly  understood.  For  this  reason  the  mathematical  modeling 
of  propellant  combustion  has  most  often  relied  on  overall  reaction 
descriptions  of  the  detailed  kinetics  in  each  macroscopic  combustion  zone. 

The  success  with  which  these  overall  reactions  can  mimic  the  effects  of 

•  I  _  O 

complex  kinetics  has  been  demonstrated  previously.  A  technique  for 
determining  the  overall  ljinetic  parameters  from  measured  temperature  profiles 
has  also  been  developed.  This  technique  utilizes  the  energy  conservation 
equation  to  compute  the  volumetric  heat  release  profile,  from  which  the 
kinetics  are  derived  in  various  ways.  Apart  from  the  spatial  derivatives  of 
the  temperature,  the  energy  equation  also  involves  the  thermal  conductivity 
and  specific  heat  of  the  local  gas  mixture.  In  the  development  of  the 
kinetics  extraction  technique  it  was  necessarily  assumed  that  these  values 
were  known.  The  present  study  was  undertaken  to  assess  the  extent  to  which 
current  uncertainties  in  these  transport  parameters  jeopardizes  burning  rate 
predictions  based  on  this  kinetics  extraction  technique. 


The  need  for  this  assessment  arose  during  a  study  of  the  feasibility  of 
molecular  beam  sampling  of  the  stable  intermediate  chemical  species  present 
during  solid  propellant  combustion.  Molecular  beam  sampling  of  radicals  had 
been  attempted  previously  but  suffered  from  low  signal  strengths.  Renewed 
interest  in  the  experiment  was  generated  by  the  desire  for  major  species 
profiles  in  combusting  nitramines  as  well  as  by  the  potential  use  in 
estimating  the  transport  properties  of  fizz  zone  gases  for  the  application 
described  above.  It  was  felt  that  such  an  assessment  would  not  only  help 
determine  if  the  molecular  beam  experiment  was  warranted  but  might  also 
produce  best  estimates  of  these  transport  parameters  for  interim  use. 

We  begin  with  a  discussion  of  experiments  in  which  gas  compositions  were 
determined  at  the  beginning  and  end  of  the  fizz  zone  of  burning  double  base 
propellants.  The  thermal  conductivity  and  specific  heat  are  then  calculated 
for  these  mixtures  and  models  for  the  temperature  dependence  of  the  properties 
are  constructed.  Concentrations  of  key  species  are  then  varied  within  limits 
thought  to  reflect  maximum  probable  experimental  uncertainties,  and  the 
mixture  properties  recomputed  for  these  extremes.  These  calculations  form  the 
basis  for  making  best  a  priori  estimates  of  the  actual  mixture  properties  and 
the  maximum  probable  errors  associated  with  these  estimates.  The  extent  to 
which  these  errors  map  into  the  derived  kinetics  parameters  and  ultimately 
into  burning  rate  predictions  is  then  determined.  (In  this  report  the  term 
"burning  rate"  will  refer  to  the  surface  regression  rate  in  either  its  linear 
form  or  mass  flux  form.) 


II.  MOLECULAR  COMPOSITION  OF  THE  PRIMARY  FLAME  ZONE 

Lengelle,  et  al.,^  have  summarized  the  current  state  of  knowledge  with 
regard  to  the  molecular  composition  of  the  gases  beyond  the  burning  surface  of 
solid  double-base  propellant.  Heat  release  in  the  gas  phase  occurs  in  two 
stages  which  here  will  be  termed  the  primary  and  secondary  flame  zones.  The 
primary  flame  is  also  known  as  the  fizz  zone  and  is  non-luminous  in  contrast 
to  the  secondary  flame.  Below  about  1  MPa,  the  propellant  burns  without  a 
secondary  flame  and  the  full  heat  of  explosion  is  consequently  not  achieved. 
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At  its  onset  the  secondary  flame  stabilizes  at  a  distance  on  the  order  of  a 
centimeter  from  the  surface.  As  the  pressure  is  increased,  the  secondary 
flame  moves  closer  to  the  surface  until  at  5-10  MPa  it  appears  to  be  in 
contact  with  the  surface.  At  these  pressures  the  combustion  zones  are  too 
compressed  to  know  if  the  two-stage  character  remains  or  if  they  coalesce  into 
a  single  zone.  The  primary  zone  appears  to  consist  of  reactions  of  NO2  with 
aldehydes,  and  the  secondary  flame  is  dominated  by  NO  reacting  with  CO  and 
H2*  At  pressures  where  the  two  stage  character  pertains,  the  preponderance  of 
evidence  suggests  that  the  burning  rate  is  controlled  by  the  primary  flame 
alone.  Accordingly,  the  focus  in  this  study  is  on  the  gases  of  the  primary 
flame . 

An  experiment  which  simulates  a  measurement  of  the  gases  emerging  from 
the  burning  surface  has  been  described  by  Lengelle,  et  al.^  ^  Under  vacuum, 
propellant  is  subjected  to  either  a  heated,  perforated  plate  or  radiant 
heating  such  that  a  surface  regression  rate  of  about  1  mm/s  is  achieved. 

These  conditions  allow  a  residence  time  in  the  condensed  phase  similar  to  that 
occuring  during  combustion  but  prevent  the  primary  gas  flame  from  developing. 
The  gases  issuing  from  the  surface  are  sampled  through  a  probe,  then  ionized 
and  analyzed  by  mass  spectrometry.  The  molecular  compositions  determined  from 
this  experiment  are  given  in  Table  1. 


Table  1. 

Surface 

Gas  Composition 

(600  K) 

ch2o 

( CHO)  2 

CO 

co2 

no2 

NO 

h2o 

h2 

xi 

X  104 

2.05e 

1.06t 

0.963t 

0.999t 

i.iot 

1 . 1 1 1 

6.95t 

Cpi 

0.383 

n 

0.780# 

e 

0.260_ 

n 

0.257n 

0.238n 

0.249n 

0.48!n 

3.48n 

Xi 

(EXP.) 

0.23 

0.11 

0.13 

0.06 

0.28 

0.11 

0.06 

0.02 

rai 

(EXP.) 

0.19 

0.18 

0.10 

0.07 

0.34 

0.09 

0.03 

0.001 

Xi 

(HIGH) 

0.46 

0.22 

0.06 

0.03 

0.12 

0.05 

0.03 

0.04 

mi 

(HIGH) 

0.37 

0.34 

0.05 

0.04 

0.15 

0.04 

0.01 

0.002 

Xi 

(LOW) 

0.12 

0.06 

0.16 

0.08 

0.35 

0.14 

0.08 

0.01 

m- 

(LOW) 

0.10 

0.10 

0.12 

0.10 

0.44 

0.11 

0.04 

0.0005 

Notes : 

X  in  cal/cm-s-K  c  in  cal/g-K 

•  12 
t  -  Touloukian,  et  al. 
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n  -  NASA-Lewis  code 

e  -  estimate  based  on  methane  (see  text) 

*  -  a  value  of  1.0  x  10  4  was  used  for  the  LOW  calculation 

#  -  a  value  of  0.30  was  used  for  the  LOW  calculation 

EXP.  -  experimental  values,'’  are  mole  fractions,  m^  are  mass  fractions 


Experiments  under  similar  conditions  have  been  conducted  by  other 
investigators  as  well.  Dauerman  and  Tajima^-^  applied  mass  spectrometry  with 
a  cone-type  sampling  port  to  the  radiative  pyrolysis  of  double-base  propellant 
at  pressures  up  to  100  torr.  Although  N0£  was  observed,  they  argued  that  NO^ 
was  the  primary  oxidizer  formed  in  the  breakup  of  the  -CONO2  groups.  Later 
work  by  Farber  and  Srivastava^®  using  mass  spectrometric  analysis  of  gases 
produced  in  heated  effusion  cells  supported  the  direct  production  of  NO2  under 
vacuum  pyrolysis  conditions. 

Lengelle,  et  al . ,  ^  have  also  extracted  samples  of  the  gases  at  the  end  of 
the  primary  flame  zone  through  small  probes  and  analyzed  them  with  mass 
spectrometry  and  gas  chromatography.  The  composition  obtained  for  propellant 
burning  at  9  atm  is  given  in  Table  2.  This  breakdown  agrees  substantially 
with  similar  measurements  made  previously  by  Heller  and  Gordon.^ 

Table  2.  Dark  Zone  Gas  Composition  (1600  K) 


n2 

NO 

CO 

co  2 

ch4 

C2H4 

h2o 

h2 

Ai 

X  10-4 

2 . 1 3 1 

2.42n 

2 . 27_ 
n 

2.46t 

6.69n 

4.80n 

C 

00 

15. 6n 

Cpi 

0.300 

n 

0.287 

n 

0.303 

n 

0.320 

n 

1.32. 

0.950n 

0.636n 

3.88n 

Xi 

(EXP.) 

0.02 

0.21 

0.38 

0.09 

0.03 

0.01 

0.20 

0.08 

mi 

(EXP.) 

0.02 

0.25 

0.42 

0.15 

0.02 

0.01 

0.32 

0.006 

Xi 

(HIGH) 

0.01 

0.11 

0.20 

0.05 

0.05 

0.02 

0.04 

0.16 

mi 

(HIGH) 

0.01 

0.16 

0.28 

0.11 

0.04 

0.03 

0.35 

0.01 

Xi 

(LOW) 

0.02 

0.26 

0.47 

0.10 

0.01 

0.005 

0.10 

0.04 

m- 

(LOW) 

0.02 

0.28 

0.47 

0.16 

0.01 

0.005 

0.06 

0.004 

Notes  : 


A  in  cal/cm-s-K  c  in  cal/g-K 

•  1 2 

t  -  Touloukian,  et  al . 
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n  -  NASA-Lewis  code 
j  -  JANAF  tables^ 


experimental  values,  are  mole  fractions,  m^  are  mass  fractions 


III.  THERMAL  CONDUCTIVITY  FOR  INDIVIDUAL  SPECIES  AND  MIXTURES 

Values  for  the  thermal  conductivity  of  individual  species  in  Tables  1  and 
2  ar|j  for  the  most  part,  obtained  from  the  data  compilation  by  Touloukian,  et 
al.,  or  the  NASA-Lewis  code.^  In  general  the  Touloukian  and  NASA-Lewis 
values  agreed  closely  with  those  computed  by  the  modified  Eucken  correlation^4 
in  which  the  viscosity  was  calculated  using  the  Lennard-Jones  potential.  In 


Xifl 


•a 


the  case  of  glyoxal,  neither  conductivity  data  nor  L-J  parameters  could  be 
found,  so  the  values  for  methane^  were  used  as  a  reasonable  upper  limit  since 
no  organic  molecule  in  the  Touloukian  compilation  has  a  higher  conductivity. 

Thermal  conductivities  for  the  multicomponent  mixtures  are  obtained  using 
a  combination  of  the  linear  and  reciprocal  mixing  rules. ^  ^ 

*mix  =  °*5  [Ei  XiXi  +  ^i  , 

where  X^  is  the  mole  fraction  of  the  species  i  and  X  ^  is  the  thermal 
conductivity  of  the  species  i.  This  empirical  equation  has  been  found^  to 
minimize  error  in  the  mixture  conductivity  since  the  linear  rule  leads  to 
overestimates  and  the  reciprocal  rule  leads  to  underestimates. 

The  pressure  dependence  of  the  thermal  conductivity  was  also 
considered.  The  estimation  procedure  of  Stiel  and  Thodos^  ^  was  used  to 
compute  the  variation  of  thermal  conductivity  with  pressure  for  the  individual 
species;  then,  the  mixture  value  was  obtained  by  the  above  empirical  rule.  At 
9  atm  the  mixture  conductivity  was  increased  over  the  1  atm  value  by  less  than 
1%.  At  50  atm  the  increase  was  less  than  5%.  Thus  only  the  1  atm  values  are 
given  in  Tables  1  and  2. 


IV.  SENSITIVITY  OF  MIXTURE  THERMAL  CONDUCTIVITIES  AND  SPECIFIC  HEATS 

TO  VARIATIONS  IN  COMPOSITION 

In  an  effort  to  determine  how  experimental  uncertainties  in  the  above 
molecular  compositions  might  affect  the  mixture  thermal  nductivity  and 
specific  heat,  we  have  arbitrarily  doubled  (and  halved)  the  mole  fractions  of 
the  species  with  the  highest  values  for  the  thermal  conductivity.  Barring 
gross  systematic  errors  in  the  experiments,  this  strategy  should  lead  to  an 
estimate  of  the  maximum  probable  error  in  the  mixture  parameters.  In  Tables  1 
and  2  the  underscored  mole  fractions  were  varied  as  just  described  and  the 
remaining  mole  fractions  recomputed  in  their  same  relative  proportions.  The 
mixture  specific  heat  is  calculated  as  the  linear  sum  of  the  raass-fraction- 
weighted  individual  specific  heats.  The  results  are  summarized  in  Table  3. 
Also  indicated  in  the  table  are  the  percent  errors  in  the  experimental  values 
should  the  extreme  values  be  correct.  It  is  clear  that  one  could  choose 
average  values  for  which  the  maximum  probable  error  would  be  less  than  25%  for 
the  thermal  conductivity  and  less  than  30%  for  the  specific  heat. 

The  variation  scheme  adopted  is  obviously  not  unique  and  is  not  even 
completely  self-consistent,  since  atomic  mass  fractions  are  not  conserved. 
However,  because  of  the  similarity  of  transport  parameter  values  among  species 
not  underlined  in  Tables  1  and  2,  attending  to  these  details  would  produce 
only  minor  corrections  to  our  analysis.  Also,  use  of  the  methane  value  for 
thermal  conductivity  in  the  case  of  glyoxal,  while  sensible  in  obtaining  an 
upper  limit  for  the  transport  parameters,  is  inappropriate  for  obtaining  the 
lower  limit.  Again,  however,  in  the  lower  limit  case,  the  mole  fraction  of 
glyoxal  is  too  small  to  influence  the  mixture  value  so  that  our  object  is  not 
jeopardized  by  this  apparently  poor  assumption.  If,  in  fact,  these  transport 
parameter  variations  did  lead  to  large  influences  on  the  burning  rate 
predictions,  then  more  care  in  their  construction  would  be  in  order.  As  will 
be  seen,  no  such  result  is  obtained. 
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Table  3.  Mixture  Properties 


Property 

Experimental 

High 

Low 

X  (cal/cm-s-K) 

@  600  K 

1.34E-4 

1.68E-4  (20%)* 

1.14E-4  (18%) 

cp  (cal/g-K'  @ 

600  K 

0.378 

0.488  (23%) 

0.278  (36%) 

X  (cal/cm-s-K) 

@  1600  K 

3.51E-4 

4.87E-4  (28%) 

2.90E-4  (21%) 

c  (cal/g-K)  @ 

1600  K 

0.395 

0.550  (28%) 

0.348  (14%) 

*  %  error  in  the  experimental  value  if  this  were  the  true  value. 


V.  EFFECT  OF  UNCERTAINTY  IN  X  AND  cp  ON  HEAT  RELEASE  PROFILES, 

OVERALL  KINETICS,  BURNING  RATE,  AND  TEMPERATURE  SENSITIVITY 

In  order  to  probe  the  effects  of  the  uncertainties  in  transport 
properties  on  the  derived  kinetics,  we  have  assembled  five  transport  models 
from  the  high  and  low  values  given  in  Table  3.  Each  model  assumes  a 
temperature  dependence  for  both  the  thermal  conductivity  and  specific  heat  of 
the  form 

f  =  aTb 

where  f  is  the  transport  property  being  considered  and  T  is  the  temperature  in 
Kelvins.  The  constants  a  and  b  for  the  five  cases  are  chosen  so  as  to 
give:  (1)  the  maximum  values  for  thermal  conductivity  over  the  range  600- 

1600  K  (labeled  MAXV  in  the  figures),  (2)  the  maximum  slope  with  respect  to 
temperature  (MAXS) ,  (3)  the  minimum  values  (MINV),  (4)  minimum  slope  (MINS), 
and  (5)  the  best  estimate  (BESTE)  of  the  values  over  the  interval.  The  best 
estimate  at  each  temperature  is  chosen  so  as  to  equalize  the  percent  error 
between  the  best  estimate  and  extreme  values  in  Table  3.  This  criterion  leads 
to  the  defining  expression  (at  each  end  of  the  temperature  interval) 

£be  ’  2  fH  <V<fH  *  fl>  ■  “> 

where  f^  and  f^  are  the  low  and  high  values  of  the  transport  parameter  f  from 
Table  3,  The  temperature  dependence  of  the  best  estimate  model  is  then 
constructed  according  to  the  power  law  function  above.  The  values  of  a  and  b 
for  all  five  cases  are  given  in  Table  4,  and  the  five  models  for  A(T)  and 
c  (T)  are  plotted  in  Figures  1  and  2,  respectively,  to  illustrate  how  they 
span  the  range  of  possibilities.  In  the  following  analysis  we  shall  adopt  the 
point  of  view  that  the  best  estimate  model  is  the  least  biased  guess  to  use  in 
the  analysis  of  fizz  zone  temperature  profiles.  We  then  suppose  that,  due  to 
errors  in  the  experimental  compositions,  the  actual  values  of  X  and  cp  m  ight 
be  given  by  one  of  the  other  four  models.  The  consequences  of  these  errors 
can  then  be  assessed  in  relation  to  calculations  based  on  the  corresponding 
transport  model. 


Table  4.  Defining  Parameters  for  A(T)  and  cp(T)  Models 


Parameter 

Min.  Slope 

Min.  Value 

Best  Eat. 

Max.  Slope 

Max.  Value 

a  ( for 

A) 

4. 7756E-6 

2.5840E-7 

2.2131E-7 

8. 7902E-9 

1 . 6245E-7 

b  ( for 

A) 

0.55659 

0.95193 

1.0037 

1.4804 

1.0851 

a  ( for 

CP) 

4.4270 

6.4258E-2 

0. 10583 

3.2469E-3 

0.22369 

b  ( for 

CP) 

-0.34472 

0.22897 

0.18876 

0.69563 

0.12194 

Note:  A  in  cal/cm-s-K,  cp  in  cal/g-K 
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Figure  1.  Models  of  Mixture  Thermal  Conductivity  in  the  Primary 
Combustion  Zone  Based  on  Maximum  Probable  Uncertainties 
in  Molecular  Compositions 

The  technique  for  deriving  overall  kinetics  from  a  combustion  wave 
temperature  profile  is  described  in  detail  in  Ref.  4.  Essentially,  it  depends 
on  the  energy  conservation  equation  simplified  by  assuming  that  the  specific 
heats  of  all  chemical  species  are  equal  and  that  thermal  diffusion  is 
negligible.  In  that  case 

(d/dx)[A (dT/dx)]  -  Mcp(dT/dx)  +  q(x)  =  0  (2) 

where  q(x)  is  the  volumetric  rate  of  heat  release  due  to  reactions  (also 
referred  to  here  as  '.he  heat  release  profile).  Equation  2  is  valid  for 


temperature  dependences  in  both  A  and  c  (though  Ref.  4  assumed  A  and  c  to  be 
constant).  First,  the  temperature  prof?le  is  least  squares  fitted  by  a^ 
suitable  but  arbitrary  analytic  function  (see  the  appendix).  The  temperature 
derivatives  in  Eq.  2  can  then  be  obtained  analytically  and  together  with 
values  for  A  and  c  one  can  calculate  the  heat  release  profile,  q(x).  The 
heat  release  profile  is  related  to  the  kinetics  by 


q  =  Qg  <PmB)  Ag  exp(-Eg/RT) 


where  Q  (cal/g)  is  the  primary  flame  zone  heat  release,  p  is  the  local  mass 


density,  mB  is  the  mass  fraction  of  our  single  hypothetical  reactant,  v  is  the 
order  of  the  overall  reaction  with  respect  to  that  reactant's  concentration, 


A  is  the  Arrhenius  preexponential  factor,  and  E  is  the  activation  energy. 

T§e  overall  kinetics  can  be  derived  from  the  heat  release  profile  as  described 
in  Ref.  4.  (The  "single  point  method"  is  used  exclusively  here.) 
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Figure  2.  Models  of  Mixture  Specific  Heat  in  the  Primary  Combustion  Zone 
Based  on  Maximum  Probable  Uncertainties  in  Molecular  Compositions 


In  this  study  the  procedure  was  applied  to  the  experimental  temperature 
:ile  published  by  Zenin1  for  a  double  base  propellant  burning  at  1.0  MPa 


profile  published  by  Zenin1  for  a  double  base  propellant  burning  at  1.0  MPa 
and  25°C.  These  data  were  chosen  because  the  temperature  profile,  surface 
temperature,  and  burning  rate  were  all  measured  in  the  same  experiment.  In 
publishing  his  temperature  profiles  Zenin  used  his  measured  value  of  the 
surface  temperature  to  locate  the  origin  of  coordinates  at  the  surface. 

Because  of  small  errors  in  digitizing  the  temperature  profiles,  the  value  of 
surface  temperature  given  by  the  analytic  fit  evaluated  at  the  surface  is 
605  K  rather  than  his  reported  value  of  about  570  K.  To  avoid  inconsistencies 
we  simply  take  the  experimental  surface  temperature  value  in  this  study  to  be 
605  K. 
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It  should  be  noted  that  the  objective  of  this  approach  is  to  obtain  an 
accurate  and  convenient  description  of  the  thermal  structure  of  the  fizz  zone 
and  not  a  detailed  understanding  of  the  chemical  kinetics.  Therefore  the 
method  should  ultimately  be  judged  by  the  fidelity  with  which  the  temperature 
and  heat  release  profiles,  computed  by  means  of  the  derived  kinetics  in  a 
single-reaction  flame  code,  can  mimic  the  original  (i.e.,  measured  T(x)  and 
derived  q(x))  profiles.  In  Ref.  4  it  was  shown  that  if  a  given  temperature 
profile  was  caused  by  a  single  irreversible  reaction,  and  if  A  and  c  for  the 
system  are  known,  then  the  kinetic  parameters  can  be  recovered  by  the  above 
procedure  with  good  accuracy.  It  was  also  demonstrated  that  overall  kinetics 
derived  from  an  actual  fizz  zone  temperature  profile,  when  used  as  input  to  a 
flame  code,  did  indeed  enable  an  accurate  reconstruction  of  both  the 
temperature  and  heat  release  profiles,  again  assuming  knowledge  of  A  and  c 
for  the  gases. 

Uncertainties  in  A  and  c  give  rise  to  uncertainties  in  the  heat  release 
profile  from  which  the  kineti§s  are  derived.  Since  the  shapes  of  the  heat 
release  profiles  based  on  the  different  transport  models  are  similar, 
differences  can  be  summarized  by  the  magnitudes  and  positions  of  the  peak  in 
q(x),  which  is  computed  from  the  fitted  temperature  profile  using  Eq.  2.  Peak 
values  and  positions  for  all  the  transport  models  are  shown  in  Figure  3.  That 
figure  indicates  that  use  of  the  best  estimate  model  results  in  the  errors  of 
the  peak  magnitude  of  probably  no  more  than  25%  but  errors  in  the  peak 
position  of  as  much  as  a  factor  of  2.  The  possible  significance  of  these 
errors  will  be  discussed  presently. 
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Figure  3.  Magnitudes  and  Positions  of  Peaks  in  the  Heat  Release  Profiles 
Determined  from  Zenin's  Fizz  Zone  Temperature  Profile  (at  10  atm) 
and  the  Various  Thermal  Transport  Models 


Table  5  shows  how  the  different  transport  models  affect  the  determination 
of  the  overall  kinetics  parameters.  Uncertainties  of  up  to  25%  arise  in  the 
activation  energy;  however,  the  significance  of  this  variation  is  best  judged 


\  s 


in  relation  to  its  effects  on  model  burning  rate  calculations  which  are 
discussed  next. 


to 

the  Temperature 

Profile 

Parameter 

Min.  Slope 

Min.  Value  Best  Est. 

Max.  Slope 

Max.  Value 

$  (cal/ctn  -s)* 

6.17  (-19%) 

4.20  (+19%) 

5.01 

4.22  (+19%) 

6.20  (-19%) 

Qg(cal/g) 

353.  (-7%) 

264.  (+25%) 

329. 

333.  (-1%) 

437.  (-25%) 

Qg(cal/g) 

171.  (-45%) 

51.  (+84%) 

94. 

52.  (+81%) 

172.  (-45%) 

E  (kcal/mole) 

o 

2.55  (+24%) 

3.32  (-5%) 

3.17 

4.12  (-23%) 

2.95  (+7%) 

V 

0.71 

0.72 

0.70 

0.73 

0.67 

Ag(cgs  units) 

4.15E3 

5.82E3 

4.79E3 

8.78E3 

3.72E3 

o 

Notes:  Based  on  the  experimental  values,  M=0.32  g/cm^-s,  T  =605  K,  and 
Tq=298  K. 

Percent  errors  are  computed  assuming  that  the  extreme  transport  models 
may  be  correct. 

*<t>  =  X  (dT/dx)  is  the  heat  feedback  flux  from  the  gas  to  the  burning 
surface. 


In  order  to  do  a  calculation  of  the  burning  rate,  a  number  of  parameters 
besides  the  gas  phase  kinetics  are  required.  We  now  discuss  briefly  how  one 
might  obtain  these  values.  Physical  and  mathematical  details  of  the 
combustion  idealization  employed  here  are  discussed  in  Refs.  19-21.  Given  a 
temperature  profile  in  the  gas  phase,  the  mass  burning  rate  M,  the  surface 
temperature  Tg,  all  determined  experimentally,  and  models  for  X  and  c  ,  the 
effective  gas  phase  kinetics  can  be  determined  as  described  above.  In 
addition,  the  gas  phase  heat  release  can  be  deduced  from  (see  eqs .  15  and  23 
of  Ref.  19,  here  generalized  for  a  temperature  dependent  specific  heat) 

T, 


Q  =  £  c  dT  +  (X/M)  (dT/dx), 

o  r  ‘ 


(3) 


where  X  is  evaluated  at  Tg  and  we  have  set  the  nascent  gas-phase  reactant  mass 
fraction  (mg®  in  Ref.  19)  to  1  in  the  absence  of  more  detailed  information. 

Tj  is  the  temperature  at  the  end  of  the  fizz  zone  and  is  here  taken  as  the 
fitted  value  at  the  last  digitized  point.  The  condensed  phase  heat  release 
may  be  obtained  from  the  energy  flux  continuity  condition  (Eq.  18  of  Ref.  19) 
at  the  gas/solid  interface: 


CP  Ts  -  cs  T0  -  (X/M)  <dT/d*), 


(4) 
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where  Tq  is  the  ambient  temperature  of  the  propellant  strand  at  the  start  of 
the  temperature  profile  experiment  and  c  is  evaluated  at  Tg .  The  condensed 
phase  specific  heat,  cg,  could  be  measured  using  standard  thermal  analysis 
methods.  Here  we  assume  a  constant  value  of  0.35  cal/g-K.  Like  the  gas  phase 
kinetics  and  Q  ,  Qg  will  be  sensitive  to  the  transport  models  as  shown  in 
Table  5.  8 

The  parameters  for  the  pyrolysis  law,  which  is  assumed  to  govern  the 
surface  regression,  may  be  obtained  from  measurements  of  the  burning  rate  and 
surface  temperature  at  a  number  of  different  pressures  and  initial 
temperatures.  This  was  done  in  Ref.  20  (NC3  data  set)  for  data  reported  by 
Zenin.  Here  we  make  use  of  the  activation  energy  so  determined  and  compute  a 
preexponential  factor  Mq,  consistent  with  the  fitted  surface  temperature  and 
measured  burning  rate  for  the  temperature  profile  used  in  this  work.  The 
pyrolysis  law  in  terms  of  the  mass  regression  rate  is 

M  =  Mq  exp  (-Eg/RTg)  .  (5) 

With  all  of  the  model  parameters  now  semi-empirical ly  determined,  one  may 
solve  the  conservation  equations  numerically  (as  was  done  in  Ref.  21)  for  the 
burning  rate,  temperature  profile,  and  heat  release  profile.  These  calculated 
quantities  may  then  be  compared  with  the  "experimental"  ones.  (Quotes  are 
used  because  the  original  heat  release  profile  requires  estimates  of  the 
transport  parameters  and  thus  is  really  a  derived  result.)  Results  of  this 
"complete  solution"  are  shown  in  Figures  4-8  and  in  Table  6  for  each  of  the 

transport  models  (all  figures  using  the  same  scale  for  easy  comparison).  In 

the  figures  the  discrete  points  are  the  digitized  values  of  Zenin' s 
experimentally  determined  temperature  profile.  The  solid  curve  through  these 
points  is  the  analytic  temperature  function  derived  by  a  least  squares  fit  to 
the  discrete  points  (see  Appendix  A).  These  two  curves  are  the  same  for  all 

five  transport  models.  The  solid  curve  with  the  maximum  is  the  heat  release 

profile  derived  from  the  analytic  temperature  profile  with  the  indicated 
transport  models.  (The  magnitude  and  position  of  the  peak  of  this  profile  are 
the  values  shown  in  Figure  3.)  The  dashed  temperature  profile  is  the  solution 
of  the  energy  equation  using  the  derived  kinetics  and  Eq.  4  as  the  solid-gas 
boundary  condition.  The  dashed  heat  release  profile  is  constructed  from  the 
overall  kinetics  and  the  dashed  temperature  profile.  For  each  transport  case 
it  is  seen  that  the  calculated  profiles  match  the  "experimental"  profiles 
quite  closely  over  the  bulk  of  the  reactive  heat  release.  The  slow  rise  of 
the  experimental  temperature  after  the  shoulder  is  characteristic  of 
relatively  slow  recombination  reactions1  and  would  require  a  separate 
effective  reaction  to  model.  Such  an  effort  is  hardly  appropriate  in  view  of 
the  small  amount  of  energy  involved. 

The  first  conclusion  to  be  drawn  from  an  examination  of  Figures  4-8  and 
Tables  5  and  6  is  that  the  ability  of  the  analysis  method  to  mimic  the 
experimental  temperature  profile,  match  the  experimental  burning  rate,  and 
give  a  self-consistent  rendering  of  the  heat  release  profile  is  remarkably 
good  for  each  model  of  the  transport  parameters.  To  some  extent  the  high 
degree  of  numerical  consistency  is  due  to  the  self-compensating  effect  of 
using  the  same  transport  model  to  solve  the  conservation  equations  as  was  used 
to  derive  the  kinetics.  Since  both  Qg  and  Qg  were  obtained  using  the 
experimental  values  of  M,  T  ,  and  Tq ,  these  conditions  constitute  something  of 
a  "calibration  point,"  i.e.,  if  the  overall  reaction  description  were 


perfectly  able  to  mimic  the  actual  thermal  structure,  one  would  expect  to  back 
calculate  the  experimental  value  of  the  burning  rate  at  the  calibration 
point.  Thus  the  very  small  errors  in  the  back  calculated  burning  rates 
provide  an  indication  of  the  basic  compatability  of  the  single  overall 
reaction  idealization  with  the  actual  fizz  zone  reaction  phenomena. 
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Figure  4.  Calculated  (Dashed  Lines)  Temperature  and  Heat  Release  Profiles 
Based  on  Effective  Kinetics  Derived  from  the  Heat  Release  Profile 
(Solid  Line)  Deduced  from  the  Analytic  Fit  (Solid  Line)  to  Zenin's 
Temperature  Profile  Data  (Points)  and  the  Minimum  Slope  (MINS) 

Transport  Model 
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Figure  5.  Calculated  (Dashed  Lines)  Temperature  and  Heat  Release  Profiles 
Based  on  Effective  Kinetics  Derived  from  the  Heat  Release  Profile 
(Solid  Line)  Deduced  from  the  Analytic  Fit  (Solid  Line)  to  Zenin's 
Temperature  Profile  Data  (Points)  and  the  Minimum  Value  (MINV) 

Transport  Model 
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Figure  6.  Calculated  (Dashed  Lines)  Temperature  and  Heat  Release  Profiles 
Based  on  Effective  Kinetics  Derived  from  the  Heat  Release  Profile 
(Solid  Line)  Deduced  from  the  Analytic  Fit  (Solid  Line)  to  Zenin's 
Temperature  Profile  Data  (Points)  and  the  Best  Estimate  (BESTE) 
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Figure  7.  Calculated  (Dashed  Lines)  Temperature  and  Heat  Release  Profiles 
Based  on  Effective  Kinetics  Derived  from  the  Heat  Release  Profile 
(Solid  Line)  Deduced  from  the  Analytic  Fit  (Solid  Line)  to  Zenin's 
Temperature  Profile  Data  (Points)  and  the  Maximum  Slope  (MAXS) 
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Figure  8.  Calculated  (Dashed  Lines)  Temperature  and  Heat  Release  Profiles 
Based  on  Effective  Kinetics  Derived  from  the  Heat  Release  Profile 
(Solid  Line)  Deduced  from  the  Analytic  Fit  (Solid  Line)  to  Zenin's 
Temperature  Profile  Data  (Points)  and  the  Maximum  Value  (MAXV) 

Transport  Model 


Secondly,  it  is  apparent  that  our  uncertainty  in  the  transport  parameters 
does  not  affect  a  model  calculation  of  the  burn  rate  significantly,  i.e.,  by 
no  more  than  two  percent.  The  heat  flux,  <J>=A(dT/dx)s,  to  the  burning  surface 
is  somewhat  more  sensitive,  in  the  range  of  20-30%.  This  latter  observation 
may  be  significant  for  use  of  the  model  under  transient  combustion 
conditions.  On  the  whole,  however,  at  least  for  the  steady-state  combustion 
case  studied,  the  effectiveness  of  the  analysis  is  not  particularly  sensitive 
to  probable  uncertainties  in  A  and  Cp. 

These  results  constitute  an  important  finding;  however,  they  only  satisfy 
a  necessary  condition  for  the  ultimate  utility  of  the  strategy.  A  more 
stringent  test  is  afforded  by  comparing  burning  rate  calculations  at  ambient 
temperatures  other  than  the  "calibration  point"  of  298  K.  Except  for  certain 
catalyzed  double-base  propellants,  it  is  generally  observed  that  burning  rates 
increase  with  increasing  ambient  temperature.  The  rate  of  increase  is 
governed  by  a  complex  interaction  of  effects  involving  both  gas-phase 
kinetics,  pyrolysis  law,  and  reaction  heats.  One  thus  expects  that  the 
different  model  parameter  sets  given  in  Table  5  will  lead  to  different  ambient 
temperature  dependencies  in  the  burning  rate.  For  these  calculations,  the 
derived  kinetics,  with  Eq.  4  as  the  boundary  condition,  are  again  used  to 
calculate  the  thermal  structure  and  mass  burning  rate.  In  this  case,  however, 
the  value  of  Qg ,  assumed  to  be  constant,  is  fixed  at  the  appropriate  value 
computed  at  the  calibration  temperature  (see  Table  5),  but  now  different 
ambient  temperatures,  Tq,  are  chosen,  and  a  corresponding  value  for  m  is 
calculated.  Figure  9  illustrates  the  different  ambient  temperature 
dependencies  for  the  transport  models  that  produce  the  largest  variation  from 
the  best  estimate  case.  As  summarized  in  Table  6,  the  effect  of  uncertainty 
in  transport  parameters  on  the  burning  rate  is  magnified  but  still  7%  or  less 
over  the  practical  temperature  range  (250-350  K) . 


»< 


AMBIENT  TEMPERATURE  /  K 


Figure  9.  Mass  Burning  Rates  Computed  as  a  Function  of  Ambient 
Temperature  Using  Effective  Kinetics  Derived  from  Zenin's 

10  atm  Data  at  300  K 


VI.  TRANSPORT  PROPERTIES  OF  THE  SECONDARY  FLAME  ZONE:  AN  ASIDE 
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In  this  paper  and  others  to  date  we  have  been  concerned  with 

describing  the  primary  flame  zone  reactions  because  the  preponderance  of 
evidence  suggests  that  the  secondary  flame  zone  does  not  affect  the  burning 
rate  at  low  pressure.  The  secondary  flame  reactions  are  of  importance, 
however,  to  interior  ballistics  since  they  release  about  half  of  the  total 
chemical  heat.  Consequently,  they  exert  a  strong  influence  on  flamespreading 
and  ignition  times.  The  method  developed  for  deriving  the  effective  kinetics 
for  the  primary  flame  is  equally  applicable  to  the  secondary  flame,  given  a 
temperature  profile  through  that  region.  The  transport  parameters  for  that 
combustion  zone  are  therefore  also  of  interest.  The  purpose  of  the  section  is 
to  establish  best  estimates  for  these  transport  parameters  for  eventual  use  in 
determining  the  effective  kinetics  of  the  secondary  flame. 


Fifer  and  Cohen^  have  reported  calculations  of  the  thermal  conductivity 
and  specific  heat  of  final  gas  phase  products  from  combustion  of  black  powder, 
HMX,  M30  propellant,  and  M9  propellant.  The  focus  of  that  work  was  on  the 
effects  of  shifts  in  the  equilibrium  species  concentrations  due  to  drops  in 
temperature  in  boundary  layers.  Results  are  presented  in  terms  of  the 
"frozen"  and  "reactive"  contributions  to  the  effective  transport  parameters. 
The  frozen  component  is  based  on  mixing  rules  applied  to  the  properties  of  the 
species  concentrations  at  equilibrium  at  some  given  temperature  and  is 
equivalent  to  the  mixture  properties  calculated  in  this  paper.  The  reactive 
component  arises  when  one  considers  steady-state  heat  flow  in  reactive  systems 
in  the  absence  of  convection  in  the  direction  of  the  temperature  gradient.  In 
terms  of  the  overall  reaction  picture,  one  has  a  single  reactant  of  mass 
fraction  m,  which  changes  according  to  the  species  continuity  equation  (Eqs.  3 
and  3,  Ref.  19), 


(6) 


(d/dx) [Dp (dm/dx) ]  +  q(x)/Qg  =  0 

Solving  for  q(x)  and  substituting  into  Eq.  2  (with  M=0),  one  gets 

( d/dx) [ (A  +  QgDp(dm/dT))(dT/dx)]  =  0  (7) 

Thus  one  recovers  the  simple  heat  conduction  equation  where  the  reacting  gas 
behaves  like  a  pure  gas  of  equivalent  thermal  conductivity  Ae££  defined  by 

A  e  f f  =  A  +  QgDp(dm/dT)  ,  (8) 

where  the  first  term  on  the  right  hand  side  is  the  frozen  component,  i.e.,  the 
true  thermal  conductivity,  and  the  second  term  the  reactive  component.  This 
equation  can  be  generalized  to  the  case  of  a  arbitrary  number  of  species 
related  by  a  full  set  of  reactions.  There  the  dm/dT  term  is  replaced  by  a  sum 
over  the  drrr/dT  for  each  species  i.  If  the  system  is  in  chemical  equilibrium, 
one  can  compute  the  reactive  component  and  treat  the  system  as  a  non-reactive 
one  with  known  effective  conductivity.  Thus  the  heat  flow  problem  can  be 
greatly  simplified  where  convection  is  not  involved  and  chemical  equilibrium 
applies.  Since  convection  is  an  important  process  in  propellant  burning,  the 
reactive  component  computed  in  Reference  23  is  not  relevant  to  our  purpose, 
though  the  frozen  component  is,  provided  the  thermal  equilibrium  assumption  is 
not  too  restrictive. 

While  the  chemical  equilibrium  assumption  obviously  is  not  valid  through 
the  secondary  flame  zone,  it  should  be  reasonable  at  the  high  temperature  end 
of  the  zone.  Presumably,  the  non-equilibrium  composition  of  Table  2  is 
appropriate  at  the  low  temperature  end  of  the  zone.  Using  Fifer  and  Cohen's 
frozen  transport  values^  at  0.95  atm  and  2957  K  for  the  double  base 
propellant,  M9,  and  our  best  estimate  values  at  1600  K,  we  expect  a  reasonable 
estimate  of  the  transport  parameters  in  the  secondary  flame  to  be  given  by 

A  =  1.49  X  10~8  t0*750  cal/cm-s-K  (9) 

cp  =  0.767  T-0,0763  cal/g-K  (10) 

for  T  in  Kelvins.  Though  the  negative  exponent  for  the  specific  heat  is  not 
typical  of  pure  gases,  in  principle  it  should  be  possible  for  shifts  in 
composition.  In  any  event  the  indicated  decrease  with  temperature  is  slight. 

VII.  CONCLUSIONS 

This  study  was  undertaken  in  order  to  assess  the  importance  of  having 
accurate  values  for  the  thermal  conductivity  and  specific  heat  for  the  gas 
mixture  in  the  primary  flame  zone  of  double  base  propellants.  These  transport 
variables  are  required  in  order  to  derive  effective  kinetics  for  the  heat 
release  in  that  zone.  Such  an  assessment  was  considered  to  be  a  prerequisite 
to  a  determination  of  the  feasibility  and  desirability  of  measuring  the  major 
species  in  the  primary  flame  by  means  of  mass  spectrometry. 


Using  published  species  concentrations  measured  for  a  double  base 
propellant,  mixture  properties  both  at  the  burning  surface  and  at  the  end  of 
the  primary  flame  zone  were  computed  and  varied  by  assuming  concentrations  of 


the  critical  species  to  be  in  error  by  a  factor  of  2.  Relative  to  a  best 
estimate  value  this  error  results  in  an  error  of  less  than  25%  for  the  thermal 
conductivity  and  less  than  30%  for  the  specific  heat.  The  relative 
insensitivity  of  the  transport  properties  to  the  species  concentrations  can  be 
understood  by  the  following  reasoning.  The  thermal  conductivity  is  roughly 
proportional  to  the  product  of  the  mean  molecular  velocity  and  the  molar 
specific  heat.  In  general,  the  larger  molecules  have  higher  molar  specific 
heats  and  lower  velocities  at  the  same  temperature.  Thus,  these  two  factors 
at  least  partially  compensate  for  each  other  as  a  function  of  molecule  size. 

The  mass  specific  heat,  being  the  ratio  of  the  molar  specific  heat  and  the 

molecular  weight  also  tends  to  vary  only  slightly  with  increasing  molecule 
size. 

The  guiding  concern  in  this  paper  has  been  the  extent  to  which  the 

probable  uncertainty  in  X  and  c  might  affect  the  accuracy  of  model 

descriptions  of  propellant  combustion.  The  sensitivity  of  model  parameters, 
derived  by  a  semi-empirical  method,  to  these  uncertainties  was  therefore  first 
determined.  The  largest  error,  about  85%,  was  in  the  condensed  phase  heat 
release  Qg ,  which  should  have  the  greatest  effect  on  the  temperature 
sensitivity  of  the  burning  rate.  The  fizz  zone  heat  release  and  activation 
energy  could  be  in  error  by  up  to  25%. 

Next  a  calculation  of  the  burning  rate  using  these  parameters  was 
performed.  It  was  found  that  the  burning  rate  computed  at  the  calibration 
temperature  was  in  error,  due  to  uncertain  transport  values,  by  less  than 
2%.  The  heat  feedback  to  the  burning  surface  could  be  inaccurate  by  up  to 
30%,  a  finding  that  may  have  significant  consequences  for  use  of  the  model  in 
transient  combustion.  Due  to  the  large  effect  of  possible  transport  errors  on 
the  condensed  phase  heat  release,  calculations  of  the  burning  rate  at  ambient 
temperatures  spanning  the  range  of  practical  interest  were  performed.  Errors 
of  up  to  7%  were  noted. 

The  conclusion  of  this  study  is,  thus,  fairly  clear.  Probable 
uncertainties  in  the  transport  properties  of  the  fizz  zone,  arising  from 
uncertainties  in  the  major  species  there,  do  not  have  a  significant  effect  on 
our  ability  to  do  accurate  model  calculations  of  the  steady-state  burning 
rate.  These  uncertainties  do,  however,  lead  to  substantial  differences  in  the 
derived  values  of  the  condensed  phase  heat  release,  which  may  be  important  for 
a  proper  treatment  of  transient  combustion,  where  accurate  descriptions  for 
the  phenomena  in  each  phase  is  more  important.  While  such  calculations  would 
be  a  worthwhile  future  task,  our  judgement  at  this  time  is  that  transport 
uncertainties  are  not  the  limiting  factor  in  the  semi-empirical  modeling 
strategy  presently  under  development,  and  therefore  is  not  in  itself  a  strong 
justification  for  doing  the  molecular  beam  sampling  experiment. 
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APPENDIX 


TEMPERATURE  PROFILE  FITTING  FUNCTION 


The  choice  of  the  functional  representation  of  the  temperature  is  based 
on  the  generally  observed  sigmoidal  shape  of  the  experimental  distributions. 
The  rationale  for  the  chosen  function  is  outlined  here. 

The  function 

S  =  [1  +  tanh(nx)] 
n 

which  in  the  limit  of  large  n  becomes  the  Heavyside  step  function, A-^  provides 
the  basic  features  of  the  required  shape.  A  number  of  parameters,  designated 
p,  through  pg,  are  included  to  give  the  function  greater  flexibility  and 
proper  scale.  Parameters  p,  and  P2  adjust  the  plateau  regions  of  the  function 
to  the  correct  initial  and  final  temperature  asymptotes,  p3  repositions  the 
function  along  the  x  axis,  p^  through  pg  form  a  polynomial  expansion  for  the 
argument  of  tanh ,  which  allows  for  asymmetry  in  the  function,  which  is  then 
raised  to  a  power,  (p^+l).  Summarizing  these  considerations: 

f ( x )  =  (p2~p1)h(u)  +  p^  (Al) 

where 

( Pq+1 ) 

h(u)  =  [(1  +  tanh  u)/2]  y  , 

„  (A2) 

2  3  4  5 

u  =  p4z  +  p5z  +  p6z  *p?z  +  p8zJ  , 


and 


z  =  x  -  p3  .  (A3) 

In  general  this  function  has  sufficient  flexibility  to  simulate  the 
experimental  temperature  profile.  Occasionally,  however,  there  is  one 
particular  region  where  the  function  does  not  adequately  reproduce  the  data. 

An  additional  expression,  having  the  form  of  a  skewed  Gaussian  function,  which 
may  be  sharply  peaked  in  one  area,  is  therefore  included.  This  expression  is 
normalized  for  greater  numerical  stability,  and  is  described  by  parameters  p^Q 
through  pj3  as  follows: 


g(v)  =  1  +  G(v)/Gn  , 

(A4) 

where 

G(v)  =  p13  exp(-a  v  +  p^v)  , 

(A3) 

a2  =  4  In  2/p122  » 

(A6) 

o 

(X 

X 

II 

> 

( A7) 

Gn  =  exp(-a2vN  +  P13vn)  , 

( A8) 

and 
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> 


(A9) 


VN  =  Pl3/2a2  ’ 

For  p^3  =  0,  i®  the  full  width  at  half  maximum  of  a  Gaussian  function 

centered  at  P^g* 

The  full  representation  of  the  temperature  function  is  then 
T(x)  =  (p2~Pi)  h(u)  g(v)  +  Pi  (A10) 

where  h(u)  is  given  by  Eqs.  (A1)-(A3)  and  g(v)  is  given  by  Eqs.  (A4)-(A9). 
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